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Abstract: For the action of an algebraic cellular automaton on a Markov 
subgroup, we show that the Cesaro mean of the iterates of a Markov measure 
converges to the Haar measure. This is proven by using the combinatorics 
of the binomial coefficients on the regenerative construction of the Markov 
measure. 



1. Introduction and main concepts 

Let (.A, +) be a finite Abelian group with unit and (^4 Z , +) be the product group 
with the componentwise addition. The shift map a : _4 Z — > A z is defined on x G A 1, 
and n £ Z by (erx)„ = x n +i. 

Let \x be a probability measure on A z and <& = J2l=e Ci<J% a nnear cellular 
automaton on _4 Z , where £,r,cg, . . . ,c r are integers such that I < r. The study 
of the convergence of the sequence ($ n (/i) : n £ N) started with the pioneering 
work of Lind in This work, devoted to the linear cellular automaton <f> = 
<j _1 + id + <t on {0, 1} Z and /i a Bernoulli measure, stated that the Cesaro mean 
of the sequence converges to the uniform Bernoulli measure. So <E> randomizes any 
Bernoulli measure. 

After such work two different strategies were developed to prove the same kind 
of results but for more general classes of linear cellular automata and other kind 
of starting measures. The combinatorial properties of the Pascal triangle and the 
technique of regeneration of measures proposed in Q] were exploited in 
to prove Lind's result in the case $ = id + a and fi is a measure with complete 
connections and summable memory decay for any finite group A. In [l^ . an 
harmonic analysis formalism was developed to prove that diffusive linear cellular 
automata randomize harmonically mixing measures (both concepts were defined in 
fljjp. Other related results were proved in [j| and 0] 

In this paper we consider the same problem for linear cellular automata defined 
in a Markov subgroup. Before state our main result we need some precise notations 
and definitions. More details and references on the subject can be found in [(tL ITi. . 
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Let © C A z be a Markov subgroup, that is, it is a topological Markov shift 
denned by a — 1 incidence matrix M and it is a subgroup of A z . For every g G A, 
let F(g) = {h G A : M gh = 1} be the set of followers of g and V(g) = {h G A : 
M.h g = 1} be the set of predecessors of g. Then T := T(fS) and V := V(0) are 
subgroups of A and F(g) and V(g) are cosets of T and V in A respectively. It 
holds that \!F\ = \V\, where | • | means the cardinality of a set. We fix a function 
/ : A — > A such that f(g) G f(g) for every g & A, and then J-(g) = f(g) + T. 

For every g e A and n > 1 define T n (g) = {h G A : M^ h > 0}. One has 
T n+1 (g) = \J heHg) T n {h) and \T n {g) \ = |.F"(0)| for every ff e A. The Markov shift 
is transitive if and only if M is irreducible, that is, there exists n > 1 such that 
!F n (0) = A. Denote by r the smallest n verifying this condition. In this case we 
have the mixing property T n {g) = A for every n > r and g G A. 

For n > and g , g n G A such that M^ o9ii > define 

£"(5o, fln) = {(5i, ■ ■ ■ , 9n-i) G ^l"- 1 : M 9i9i+1 = 1, i G {0, . . . , n - 1}}. 

Then |£ n (.g ,.g„)| = |£ n (0,0)|. Therefore, in the transitive case, \<L n {g,h)\ = |£ n (0,0)| 
for every g,h <E A and n > r. 

The Haar measure of (<S, +) is denoted by v. It is the Markov measure defined 
by the stochastic matrix L = {L g h ■ g,h G A) with L g h = \J-\~ 1 if h G ^(g) and 
= otherwise, and the L— stationary vector p = (p(g) = |^4| _1 : g G A). The 
Haar measure is the maximal entropy measure for the Markov shift (©, a). 

It is useful to introduce a notation for the finite-dimensional distributions asso- 
ciated to v. Thus, for I > 1, is the finite-dimensional distribution of v in A e , 
and it is concentrated in the subset &£ — {(go, ■ ■ ■ , ge-i) G ^ : 5,+i G J-(gi), Vz G 
{0,...,£-2}}. Hence, 



l.J. Main Result 

Let jU be a Markov probability measure in ^l z that is defined by a stochastic matrix 
P and by some probability vector ir invariant for P (it ■ P = 7r). The measure ^ is 
said to be compatible with M if P g h > if and only if M. g h > 0. This property is 
equivalent to the fact that the support of a is ©. 

If x G A Z and m, n G Z are two integers with m < n, define 



Consider the endomorphism $ = id + a on _4 Z , that is ( < E>x)„ = + x n+ \ for 
x G A 7 * and n G Z. Our main result is the following one. 

Theorem 1.1. Assume a is a Markov probability measure in A z compatible with 
M. Furthermore, assume the Abelian group A is p s -torsion for some prime number 
p and s > 1 and that M is irreducible. Then the Cesdro mean of pi, under the action 
o/<f>, converges to the Haar measure v. That is, for any m G N and g G <& m : 



if 7£ 





(1.1) 
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The elements of our proof includes a regenerative construction of the Markov 
measure and the combinatorics of the binomial coefficients associated with the 
iterates of the cellular automaton. In 0] a more general result is shown: for Markov 
fields verifying a 'filling property' it is shown that the attractive property of the 
Theorem holds. In the one-dimensional case this 'filling property' is always satisfied, 
therefore the result follows. 

2. Construction of a Markov measure and renewal properties 

We use the Athreya-Ney pj representation of Markov chains. It says that it is 
possible to enlarge the probability space we are considering in order to include a 
family of integer random times. 

2. 1 . Construction of a Markov measure 

Let © C A z be a Markov subgroup with incidence matrix M. Let \x be a Markov 
probability measure in A z compatible with M defined by a pair (tt, P) as in sub- 
section 11.11 We describe a procedure to construct the restriction of /i to _A N . In 
this purpose it is useful to introduce the following notation: for g G A we put 
\i g the measure induced on _4 N by /i conditioned to the event = g}; then 

X) s e^4 7r (5)/ i g coincides with the restriction of /i to _4 N . 
Let a > be a strictly positive number such that 

a < mm{P gh :g,heA, M gh > 0}. 

We consider a probability space (f2, B, P) and three independent processes of i.i.d. 
random variables defined on this space U = (U n : n G N), W = (W n ■ n e N) and 
V = (V n : n <E N) whose marginal distributions are as follows: U n is Bernoulli (a), 
that is P([/ n = 1) = a = 1 — P(U n = 0); W n is uniformly distributed in J 7 , so 
P(W / T! = g) = \J~\~ 1 for g G T; and V n is uniformly distributed in the unit interval 
[0,1]. 

Let us construct a sequence (x n : n e N) G .4 N as a deterministic function of 
(U, W, V). For any g,h G A define 

Pgh-al^Mgh 
I4gh — : • 

1 - a 

Thus, Q = (Q g /i : g,h E A) is a. stochastic matrix compatible with M. For each 
g G A fix {Q g h Q [0, 1] : h G ^F{g)} a measurable partition of the interval [0, 1] such 
that the Lebesgue measure of Q g h is Q g h- For all g G A, u G {0, 1}, w G T and 
v G [0, 1] define 

H(g,u,w,v) = u(f(g)+w S J + (1-u) ^ h\~^{v). 

her(g) 

Now, for any G A and for n > we put 

x n = H(x n -i,U n ,W n ,V n ). 

It is clear that the distribution of the sequence x = (x n : n G N) is n x _ 1 - If £-i G A 
is a random variable with distribution n, then the distribution of x is \x. 
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2.2. The associated renewal process 

For any s, t G N with s < t define 

( U* = 1 ) ^=> ( U k = 1, for all k e {«, t} ) . 

For every m > 1 define a renewal process (Ti m) : n G N) given by: 
T (m) = 0, 7\ (m) = min{z > T (m) : Uj +m = 1}, 

and 

= min{i > T^ m \ + m : U* +m = 1} for n > 2. 

It is clear that (T^\—T n m ^ : n > 1) is a family of i.i.d. random variables. Also, from 
our computations below, it follows that the distribution of has a geometrical 
tail; that is, there exists (3 := (3{m) G [0, 1) such that P(T 1 (m) > t) < (? for any 
t>0. 

Let N^™ 1 ) be the renewal process induced by [T n m ^ : n G N). That is, for any 
iCN, 

N( m )(^) = { n e ^ : for some feN, T/ m) = n}. 

Let n G A with n > 0. One has that ?7„_i = and U™ +m = 1 implies n G 
N(™)(A). Then P(n G N( m )(A)) > a m+1 (l - a) := S > 0. Clearly if G A then 
P(0 G N( m )(A)) = 1>5. 

For A a finite subset of N and m G N we say j4 is m-separated if, for any a,6e 4 
with a^i), |a — 6| > m + 1. One gets, 

( A is m-separated ) => ( P (N( m )(A) ^0) > 1-(1-(5)I A I ) . (2.1) 
Let A( ro ) be the largest m-separated subset of A; then |yl( m )| > \A\/(m+l). Thus, 
P ( N^(A) = ) <p(N( m )(A( m )) = 0) < ((l-^ 1 ^"^) 1 



Hence, the distribution of has a geometric tail. 

3. Convergence of the Cesaro mean 
3.1. Independence lemmas 

Assume is transitive and recall r is the smallest integer verifying !F r (0) = A. 
Also recall that 7£ = |.4| _1 |.F| _ ^ _1 ). Let m < n in Z. If x is a random variable in 
A z with distribution /j, define J 7 ^ to be the sigma-algebra generated by x^. 

Lemma 3.1. Let k > r and m > 0. Then the random variable x^ +m conditioned 
to (u£+" 1 = l,^J- r - 1 ) is v {m+ ^- distributed. That is, for any g G <5 m+ i, 

r ^ x fe g U fc-r — -^0 y— 7m+l- 

,4Zso, the variable x k + m conditioned to (U*+™ +r = 1, J 7 ^ 1 V Fg +m+r+1 ) is 
i/( m+1 ) — distributed for any n> k + m + r + 1. That is, for any g G <3 TO +i, 



( k+m _ 
[ k - 



•w-rk+m+r i sr-k — r—1 w <rn \ 
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Proof. Let g = (go,...,g m ) € © m +i and put n = k — r. For any fixed h 

(ho, . . . , hn-l) G ®n, 



x£ +m = g x"- 1 = h and U£ +m = 



(x£ _1 = z and x fe+m = g 



x"" 1 = h and Ut +m 



»€e r + 1 ()in-iiflo) 

53 = 1^+1(0,0)1 • |jT|( r+m+1 ) 

aeC r+1 ()in-i.ffo) 

|T|r+l i 
r I I -r-l-fr+m+l) 1 



K I 1/11 I cm I 7m+l 



|^| ' ' \A\-\J* 

This proves the first part, the second one is entirely analogous. □ 

Now assume that A is p s -torsion for some prime number p and some s > 1, with 
s being the smallest number verifying this property. That is: mg = for every g G A 
and to G p s Z, and for every to < p s there exists some g G A such that mg =/= g. 
Observe that for every c G Z pS relatively prime to p, there exists a multiplicative 
inverse c _1 G Z pS , such that cc _1 = f mod (p s ) and c _1 is also relatively prime 
to p. Thus, for any g G A, ccr x g = g. Moreover, (eg = h) (g = c h). 

3.2. The transformation 

Recall a is the shift map in A z and $ = id + a is an endomorphism of A z . Fix 
x G -4 Z . Then for all n > and igZ, one has 

(#"x)i = (3.1) 

For every to G N denote by to/ s ) its equivalent class mod (p s ) in Z p «. Hence, for 
all n > and i 6 Z, equality l|3.1|l can be written 

n , v (s) 

(<J>"x) i = £ x i+fe . (3.2) 



Let to,£ > 0. For n > and < k < n we say fc is (to, £) — isolated in n if (?) is 
relatively prime to p, while, if fc' G {fc — to, fc+£} with k' k then (1/)^ = 0. 
Here the convention is (3) — whenever fc' < or fc' > n. 

Lemma 3.2. Letm G N andn > 2r+2m+l. J/0 < fc < n is (r+m,r+m) -isolated 
in n, then for every i G Z and g G (5 m +i, one nas 



($n x)r m = g U^^r" 1 = 1J = 7m+l. 

Proof. Since /x is cr-invariant then $ n (^) is <r- invariant too; hence it suffices to prove 
the result for i — 0. In other words, it suffices to show that 

p f($"x)» = g u*t;+™ = i) = 7m+1 . 
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(s) k — r — m—1 / \ (s) n , \ (s) 



*i = ELj w ft e y w+ e y w> 

(3.3) 

where (*) is because k is (r + m, r + m)-isolated. 

Thus, ($ n x)j = X,j + Yj. If X = pf , . . . , X m ) and Y = (Y , . . . , Y m ), then 
($™x)™ = X + Y. One gets, 



($™x)™ = g u£+; +m = i)=p(x + Y = g 



P (y = g - h 



^k+r+m = 1 



X = h, \j k + r + m 



k—r 



x P X = h 



U 



(3.4) 



Let c = (?)^ j then c is relatively prime to p, and Y = c ■ xj; +m . Thus, if c 1 

fcH 
fc 



is the (mod p s ) inverse of c, then x^ +m = c 1 • Y. Thus, for any fixed h £ A m+1 , 



Y = g-h 



X = h and U<£+ r +™ = l 

K—r 



(3.5) 



(t) 



x^- = c- 1 (g-h) 



X = h and U 



where (f) is because equation (|3.3|) implies that X is a function only of x r 1 and 
and it allows to apply Lemma Substituting I J3.5J I into (|3.4(1 yields 



(($"x)y 



Therefore the result follows. 



!] = E 7m+i-P(X = h 

= 7m+l- 



u 



/c+r+m 
k—r 



□ 



3. 3. Elementary facts on the Pascal triangle 

Now we use the following result on the Pascal triangle. Let n = J2jen n jP'' ^ e ^ ne 
decomposition of n in base p, so nj £ {0, . . . ,p ~ 1} for every j £ N. For < k < n 
consider the decompositions k = YliEN^jP 1 an d n — i 



Kummer's Theorem on binomial coefficients, whose proof can be found in [4|, states 
the following result. 

Lemma 3.3. The biggest integer i such that p e divides u) is the number of carries 
needed to sum k and n — k in base p. 

We introduce the following notation. For n £ N and i > we put 

Ji{n) = {j>i : nj^O}, &(n) = \Ji(n)\. 
For a real number c denote by [c\ the integer part of c and for a > define 
p a m = { p a ■ n : n £ N} = {n e N : = for all < j < a}. 



The 
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Lemma 3.4. Let to > 1 and a > 2s + l be such thatp^ a / 2 ^ > m. For every n G p a N 



and i > a, 



{0 < k < n : k is (m, to) -isolated in n} 



> 2^+*(") - 1. 



Proof. Fix a nonempty subset J C J i+a (n), and define < k < n by kj = rij if 
j ^ J and fcj = rij — 1 for j S J. Therefore n — k verifies (n — fc) j = 1 for j G J and 
(n — fc)j = for j ^ J. Then there is no carry in the sum of k and n — fc, so Lemma 

13.31 savs that (?) is relatively prime to p. It remains to show that (?,) = for all 
k' € {k - to, . . . , k + to} \ {fc}. 

Case 1: Let 1 < v < to and fc' = k — v. Then the p-ary decomposition of k' has 
some nonzero elements in coordinates between and a — 1 (because p a divides k, 
but does not divide v)\ moreover, it has at least b — \_a/2\ zeros in {0, . . . , a — 1}. 
However, the p-ary decomposition of n has no nonzero elements in {0, . . . , a — 1}, 
so there must be at least b carries in the addition: (n — k') + k' — n. Thus, Lemma 
EDsays that p b divides (™). But b = [a/2\ > [(2s + 1)/2J > s, so we conclude that 

G?) w 

Case 2: Let k! = k + v with 1 < v < to. Then k'^ > nj for every j < a + i and for 
some j' < a we have k'y > ny . Hence, the sum in base p of k' and n — k' will have 

at least a carries, so Lemma l3~3l savs that p a divides (2) and finally (2) =0. □ 

Lemma 3.5. Let a > 0. Then, the set 



M = 



[nep a N : e a+L , logp(n) j(n) > Jlog» 



is of density 1 in p a N. That is, 

hm l^on{Q,-^-i}| = L 

jv^oo |p a Nn{0, ...,N-l}\ 
Proof Let M = [n e N : ^ log »j (») > §(1°S» + «)}• Then 

Mn{Q,...,iV- 1} =1. 



lim — 

iV->oo iV 



To see this, let n s N be a 'generic' large integer; then the Law of Large Numbers 
says that only about ~ of the p-ary digits of n are zero; hence are nonzero. 
Since there are ^log p (n) digits in the range log p (n), log p (n)], we conclude, with 
asymptotic probability 1, that at least log p (n) digits in [| log p (n), log p (n)] are 
nonzero; hence logp(n) j (n) > ^-\og p {n) > g(log p (rc) + a) (assuming p > 2 
and log p (n) > 4a). 

Now define bijection i\> : N — > p a N by ^(rz.) = p a rt; then ^(.M) 3 M - The 
lemma follows. □ 



3.4. Proof of Theorem [7771 

Proof. Fix to > 1. If TV C N, we say that the TO-dimensional marginal of the Cesaro 
mean converges along J\f if for any g£(B m 



lim - — - r — — 

iv^oo \Afn{0,...,N-l}\ 

1 1 neMn{Q,...,N-i} 



^ ( rxir 1 = g ) = Kg). (3.6) 
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Let m' = m + r and consider a > 2s + 1 with pL a / 2 J > m > ( as i n Lemma f3.4|) . Let 
A^o be as in Lemma I3~5l We claim that the m-dimensional marginal of the Cesaro 
mean converges along Ada- 

Let n G Aio be enough large such that i = [|log p (n)J > a. Define A = 
{0 < k < n : k is (m',m') -isolated in n }. Therefore, by Lemma 13.41 and the defi- 
nition of A^o, 

|^| > 2 ^+i(«) - 1 > 2S log p ( ") -1 = 2* l ° 6 ^ n) - 1 = n 1/c - 1, (3.7) 
where C = 5 log 2 (p).Thus, 

(i (3k G A with U£+™' = l) > /i (n( 2 " 1 '' (A - m') ^ 0) 

> > 1-(1- J)" 1 ' " 1 , 

where "l|3.7jl" is by equation Ij3.7jl and " l|2.1jl " is by equation l|2.1jl (since A is 
(2m')-separated). 

Finally, LemmaP implies |/i(($ n x)| +m_1 = g)— 7m | < (l-S) nl/C -\ and thus, 

lim^oo.neMo |/x(($™x)^ +m_1 = g) - 7m| = lim, woo (l - 5)" I/C_1 = 0, as desired. 

Lemma 13.51 then implies that the m-dimensional marginal of the Cesaro mean 
converges along p a N. Now, since v is invariant for powers of <!>, we find that for any 
< j < p a , the m-dimensional marginal of the Cesaro mean also converges along 
M.j = {n + j;n G Aio}- Therefore follows from the fact that 

1 N ^ 

= E |*n{0, ■!.,*-!>! E /((W= g ) 

0<?<P a 1 L J 1 neMjn{0,...,N~l} 

= Kg) □ 
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